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We employ the recently introduced conformal iterative 
construction of Diffusion Limited Aggregates (DLA) to study 
the multifractal properties of the harmonic measure. The 
support of the harmonic measure is obtained from a dynam- 
ical process which is complementary to the iterative cluster 
growth. We use this method to establish the existence of a se- 
ries of random scaling functions that yield, via the thermody- 
namic formalism of multifractals, the generalized dimensions 
D q of DLA for q > 1. The scaling function is determined just 
by the last stages of the iterative growth process which are 
relevant to the complementary dynamics. Using the scaling 
relation D3 = Do/2 we estimate the fractal dimension of DLA 
to be D = 1.69 ±0.03. 

The diffusion limited aggregation (DLA) model was in- 
troduced in 1981 by T. Witten and L. Sander [§. DLA 
has attracted enormous amount of research as an elegant 
example of the dynamical creation of a non-trivial fractal 
set, and as a model that underlies many pattern form- 
ing processes including dielectric breakdown 0j , two-fluid 
flow H, and electro-chemical deposition In spite of 
the significant amount of effort to understand the frac- 
tal and the multifractal properties of DLA, there exists 
to date no accepted calculation of these properties from 
first principles. In this Letter we present a theory that 
is based on an iterative conformal construction of DLA 
||fr| that culminates with the construction of a series of 
(random) scaling functions which allow, via the thermo- 
dynamic formalism of multifractals , a convergent 
calculation of the multifractal properties of DLA. 

Originally the DLA model was introduced as the out- 
come of N random walks. Fixing one particle at the cen- 
ter of coordinates in e?-dimensions, and releasing random 
walkers from infinity one at a time, one allows them to 
walk around until they hit any particle belonging to the 
cluster. Upon hitting they are attached to the growing 
cluster. We are interested here in d = 2 for which nu- 
merical simulations indicated that for N — > 00 the clus- 
ter attains a fractal dimension of about 1.71 ||. The ap- 
proach used here is different: we employ an iterative con- 
formal construction of DLA that was recently proposed 
by Hastings and Levitov ||. The basic idea is to follow 
the evolution of the conformal mapping $W ( w ) which 
maps the exterior of the unit circle in the mathematical 
w-plane onto the complement of the cluster of n parti- 
cles in the physical z-plane. The unit circle is mapped 
to the boundary of the cluster which is parametrized by 
the arc length s, z(s) — & n \e %9 ). This map & n \w) is 
made from compositions of elementary maps 4>\ 8, 



$ (n) (w) = $ 



- (bC™- 1 ) 



(1) 



where the elementary map (j)x,e transforms the unit cir- 
cle to a circle with a "bump" of linear size VA around 
the point w = e 10 . Accordingly the map ^ n \w) adds 
on a new bump to the image of the unit circle under 
$(«-!) ( UI ), xhe bumps in the physical z-plane simulate 
the accreted random walkers in the original formulation. 
The main idea in this construction is to choose the po- 
sitions of the bumps n and their sizes v'An such as to 
achieve accretion of fixed linear size bumps on the bound- 
ary of the growing cluster according to the harmonic mea- 
sure P(s)ds. The latter is the probability that a random 
walker would hit an infinitesimal arc ds centered at the 
point z(s). This is done as follows. 

The probability density P(s) is given by the in- 
verse of the derivative of the conformal map P(s) — 

$' (n) (e ze ) . Using the obvious fact that ds = 

\^ n \e ie )\d9 we see that P(s)ds = d0. i.e. the harmonic 
measure is uniform on the unit circle. Thus choosing 
random positions n , and A„ in Eq.(|l|) according to 



|$(n-l)'( e ie„)|5 



(2) 



we accrete fixed size bumps in the physical plane accord- 
ing to the harmonic measure. Finally the elementary 
map (j)\,e is chosen as 



4>\,o{w) = w 1 



(1 + A) 
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4>\,9(w) = e <j)\fi(e w) . 



(3) 
(4) 



The parameter a is confined in the range < a < 1, 
determining the shape of the bump. In this Letter we 
employ a = 2/3 which is consistent with semicircular 
bumps. The qualitative properties of this mapping that 
enter prominently our analysis are the following: 



<P\,e(w) fa (1 + X) a w for 



1 



Va 



> const , 



(5) 



where const here is a tolerance-dependent but A- 
independent constant of the order of unity. 

\<f>\,e(w)\ m (l + VAJH for argw e [0 ± y/\] , (6) 
\Mw)\ » \w\ for argu; g [6 ± VA] . (7) 

Eq. (^) means that points that do not belong to the bump 
are only reparametrized along the circle. 
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The recursive dynamics can be represented as itera- 
tions of the map 4>\ n ,e n {w), 



(8) 



The final calculation of A„ can be done many times (for 
a given history 61,62,..., #n-i) according to the uniform 
measure on the unit circle, yielding a distribution of val- 
ues of A„, which according to Eq.(||) is the distribution 
of the density of the harmonic measure. The moments of 
this distribution are connected to the generalized dimen- 
sions of the harmonic measure. The latter are defined as 
follows H : consider a partition of the cluster boundary 
into balls of diameter £i and measure pi. The generalized 
dimensions D q are determined by the equation 



M 



Pi 



lim / ..„ 



--1 H 



(9) 



The calculation of D q from the statistics of X n was dis- 



cussed in detail in Ref. jjj with the exact result 



-2qD 2q+1 /D 



(10) 



It is well known ]lQ| that the fractal dimension is Dq in 
this language, the information dimension (known also as 
the dimension of the harmonic measure) is D\ = 1 JlT| , 
and in general D q > D q i for any q' > q. It was shown in 
|L2| and in @ that D 3 = D /2. This last result means 
of course that on the average A„ decreases like n -1 . The 
exact result found in is 



An — 



aDn 



(11) 



Consider the support of the harmonic measure. We use 
the uniformity of the measure on the unit circle to form, 
for a given cluster of n bumps, an equi-measure partition 
of M balls in the physical space, by selecting points z k 

z k = $(»)( e ^) = Al(fll o A2 ^ 2 o . . . o 4>x n ,eAe l(k ) , 

(12) 

where the points k = 1, 2, . . . , M are uniformly spaced 
on the unit circle. We introduce the "complementary 
dynamics" to the cluster growth by 



■ </'A„,e„(e lCfc ) 



(13) 
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FIG. 1. The evolution of the set {zj^}. (a) j=ll and 1, 
(b) j = 10 4 , 6000 and 2000. Note that the straight lines are 
not parts of the set, they simply connect points on the set. 

This observation is easy to understand. As noted in 
Eqs.(Uj^) the elementary map <p\.e(w) distorts circles of 
radius close to unity, but acts as a uniform multiplica- 
tion for points with larger absolute values. For j = n we 
always start the complementary dynamics on the unit cir- 
cle, and (|)-0) are applicable. We estimate now how many 
iterations of the complementary dynamics are necessary 
before (^|) becomes applicable. Consider an arbitrary 
value of Cfc . The first iteration of the complementary dy- 
namics z ni k have \z n ,k\ ~ 1 + iff Cfc S [6>„ ± y/K], an 
event of probability n n = \/\^/2ir. Otherwise \z n .k\ = 1. 
If the dynamics failed to increase |<z n ,fe|, it can do it with 
l z n-i,fc| after two steps, with probability (1 — k„)r„_i. 
The average number of steps A needed to grow out with 
certainty is therefore 
00 

A=y]k(l- K n )(l - K n -i) ... (1 - K n _ k+2 )K n -k+l 



k=0 



(14) 



Averaging over the history Ai,...,A n , approximating 
(A„_fcA„_ fc /) - (A„_ fe )(A„_i/) (which was justified for 
k =/= k', k, k! < n in M), we find 



and zi^ = Zk- The complementary dynamics creates 
the points z k from the seeds e l( ' k . In Fig.l we display J2, 

the typical evolution of the set {zj^} for different values (^) = ^(1 — ••■(!— (nn-k+2)){K n -k+i) (15) 
of j. The striking observation is that the equi-measure fe=0 

partition is fixed in shape very early (large j) in the com- We assume and show self-consistently that the sum is 
plementary dynamics. Later stages only serve to inflate dominated by k < n. In that case 
the fixed shape, with the last steps being most prominent 
in determining the final radius of this set which of course 
is the radius of the cluster. 
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,D 2 /D 



, (16) 
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where ([H]) has been used. The variance is estimated 
similarly, and is of the same order. Since 0.5 = D3/D < 
D2/D < Di/D — 0.58 we see that for large n we need 
relatively few iterations of the complementary dynamics 
to increase the radius of the set to a value after 

which the overwhelming majority of the iterations serve 
simply to inflate the radius by factors of (1+Xj) a . In light 
of Eq.([ll|) we understand the phenomenon exhibited in 
Fig. 1 , that the last steps of the complementary dynamics 
have the largest effect on the radius of the set {zj t k}- We 
note in passing that these comments also explain the find- 
ing of |j]|7| that the first Laurent coefficient of 3>( n \w), 
denoted there as F± is a measure of the radius of the 

cluster. The exact result is that Fi = n"=i(l + Ai) a . 
Finally note that the set {zj t k} for fc — * 00 generates 
the support of the harmonic measure. It will reveal the 
generalized dimensions D q for q>l only. 
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FIG. 2. (a) The random scaling function P2(cr) obtained 
from two consecutive partitions of 16 and 32 equi-measure 
balls. 550 clusters of 10000 bumps were used, (b) P2(cr) 
computed from 32 and 64 balls, with 550 clusters of 10000 
bumps, (c) Pi{cr) computed from 32 and 64 balls, with 44 
clusters of 100 000 bumps. 

At this point we can introduce a "binary" scaling func- 
tion. Consider two values of M, M\ = 2 m and Mi = 
2 rn+1 . Consider the two associated sets {z i } 1 ^} 1 , {zk}^^ 
such that Zi — zn. For every difference Zi — ii-i in 
the coarser resolution Mi (denoted as the ith "mother"), 
there are two differences zn — Z2i-i and Z2i-\ — ^2i-2 in 
the finer resolution, which are denoted as "daughters" . 
The binary scaling function at every resolution has M 2 
values which are obtained as the ratio of daughters to 
mothers, 

<J2l = (\%2i - Z 2 i-l\)/(\Zi - Zi-l|) 
<?2i-l = [\Z21~1 ~ Z 2 i-2\)/(\Zi - ■ (17) 



Conjecture: For 1 <C M2 <C n the binary scaling 
function ( ji~7j ) converges in distribution to a universal 
function independent of M2 and the history 6\,...6 n . 
We denote the distribution as ^(c). 

In Fig. 2 we show the numerical evidence for this con- 
jecture. We employ clusters with n — 10 4 and n = 10 5 , 
and various values of M2. Note that M2 <C n; increas- 
ing M2 beyond, say, 128 for n = 10 5 results in exposing 
the ultra-violet cutoff of the smooth bumps, yielding a 
spurious peak at a = 1/2. The harmonic measure is ex- 
tremely concentrated near the tips of the cluster, and any 
attempt to resolve the fjords leads to oversampling of the 
smooth bumps on the tips. Note that the scaling func- 
tion is defined as a ratio, and the discussion above implies 
that it is sensitive to only n D2 / D last growth steps. The 
majority of the iterations of the fundamental map <p\ t g 
are irrelevant, as they cancel in the ratio. 

Ordering the 2" diameters Zj — Zj-i in a binary basis 
£(e m . . . ci), redefining the ratios a(e m . . . ei) accordingly, 
and using in Eq.(||) the fact that pi — 2~ m for the Mi 
balls, we derive the equation Q 

<r- T{g) (tm+i...ei)r T M(e m ...e 1 ) 

e m +i...ei 

= 2 « r T W(e m ... ei ) , (18) 

e m ...ei 

where r(q) = (q — l)D q . Iterating, we find 

S m+1 = ^ T(q) • ■ ■ ti)°~ T(q) (£m • • • ei) • ■ ■ <J~ T ^ (ei) 

£ TO+ i...ei 

= 2 {m+l ^ q . (19) 

To compute Di we notice that r(l) = 0. For r(q) — > + 
all the realizations of the products of random numbers 
a~ T are comparable. Since the most probable product 
for m — * 00 is exp[(m + l)lner], we can estimate the sum 
in |TJ) as 2( m+1 > exp[(m + 1)W]. Substituting in @ 
and taking the mth root and the log to base 2 gives 

hm r(q) = . (20) 

9^1+ log 2 a 

Using the scaling function P2 (a) wc compute log 2 a ~ 
— 1 ± 0.01, yielding the expected result D\ = 1. This is 
the first nontrivial calculation of a generalized dimension 
in this approach. 

For q > 1 the approach using the most probable prod- 
uct of random numbers is not applicable. We have 2 m 
realizations of products of m random numbers, and rare 
events are relevant. Moreover, one should notice that 
even though every factor of o~ T is independent, the prod- 
ucts in the sum are not, as they have common factors. 
One can use the fact that the random products are orga- 
nized on a binary tree to write an exact recursion relation 

S m+1 = *(0)-tM(P)S<» + a(l)-^5« , (21) 
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where Sm are two independent realizations of the sum of 
2 m products, each of m random variables. We failed to 
find the exact asymptotics of S m which is necessary for 
computing r(q). We recognize however that our binary 
partition is in fact arbitrary, and instead we can perform 
refinements into k daughters at each step. Accordingly 
we will have fc scaling ratios cr(e) for each mother, with 
e now taking on k values k = 0, 1, . . . k — 1. The random 
scaling function is now denoted as Pfe(cr); its existence as 
a universal function for k = 2 m emanates in an obvious 
fashion from the existence of Pii^), and its existence for 
any k can be demonstrated independently as done above 
for P 2 (a). The important point is that the asymptotics 
of Eq.(|19|) is computable in the limit k — > oo: 



The best fit predicts 2D 3 = 2D 3 (k -> oo) = 1.69. Exam- 
ining different partitions of the unit circle and different 
integration schemes we concluded that we can bound the 
errors around D = 2D 3 = 1.69 ± 0.03. 

The excellent result of this calculation leaves for fu- 
ture research the analytic determination of the random 
scaling functions Pfc(cr). If these could be written down 
from first principles the problem of DLA will be settled. 
In particular only analytic Pk(cr) will allow reliable cal- 
culations of D q for q — > oo due to the sensitivity to the 
small a tail where the statistics is low. 
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£1=0 

In the limit of k — > oo this equation reads 



S, 



m+l 



ka(e\ 



r {q) <? 



-T(«) 



(22) 



(23) 



is random, 



where we have used the fact that cr(ei) 
independent of the consecutive factors a{e\,e2) T ■■■ 
which consist . Asymptotically 5 m +i — > 

[ka^ T ^] m+l , and substituting in ( |l9| ) we compute r(q) 
from 



P k {a)o-- f ^do- = k q - 



lim f(q) = T(q) . (24) 



FIG. 3. 2Dz(k) as a function of 1/fc. We are interested 
in Do = 2Di(k — > oo). The best fit (dashed line) predicts 
D =1.69. 

In Fig. 3 we show the solution 2D 3 — f(3) from j2^ ) 
as a function of 1/fc. It should be noted that when fc in- 
creases the determination of f(q) becomes easier since the 
LHS of ( p4| ) changes from very small to very large value 
over a narrow range of value of f . This is compatible with 
the thermodynamic formalism of multifractals p3[ . We 
expect the limit as 1/fc — > to be regular, and we fit a 
3rd order polynomial to the points as shown in the figure. 
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